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On the three-dimensional finite Larmor radius 
approximation: the case of electrons in a fixed background 

of ions 

Daniel Han-Kwan* 
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^ ' Abstract 

o ■ 

This paper is concerned with the analysis of a mathematical model arising in 
plasma physics, more specifically in fusion research. It directly follows |17j . where 
, the tri-dimensional analysis of a Vlasov-Poisson equation with finite Larmor radius 

scaling was led, corresponding to the case of ions with massless electrons whose density 
follows a linearized Maxwell-Boltzmann law. We now consider the case of electrons in 
Ph ' a background of fixed ions, which was only sketched in [T7]- Unfortunately, there is 

. evidence that the formal limit is false in general. Nevertheless, we formally derive a 

fluid system for particular monokinetic data. We prove the local in time existence of 
analytic solutions and rigorously study the limit (when the Debye length vanishes) to 
a new anisotropic fluid system. This is achieved thanks to Cauchy-Kovalevskaya type 
techniques, as introduced by Caflisch [7] and Grenier [T^. We finally show that this 
approach fails in Sobolev regularity, due to multi-fluid instabilities. 

Keywords: Finite Larmor Radius Approximation - Anisotropic quasineutral limit 
^ . - Anisotropic hydrodynamic systems - Cauchy-Kovalevskaya theorem - Ill-posedness in 

Q ! Sobolev spaces. 

1 Introduction 
1.1 Presentation of the problem 

The main goal of this paper is to derive some fluid model in order to understand the 
behaviour of a quasineutral gas of electrons in a neutralizing background of fixed ions and 
' submitted to a strong magnetic field. For simplicity, we consider that the magnetic field 

has fixed direction and intensity The density of the electrons is governed by the classical 
Vlasov-Poisson equation. We first introduce some notations: 

Notations. Let (ei,e2,e||) be a fixed orthonormal basis ofM.'^. 

• The subscript _L stands for the orthogonal projection on the plane (61,62), while the 
subscript \\ stands for the projection on ey . 

• For any vector X = (Xi, X2, Xy), we define X-^ as the vector {Xy, —Xx, 0) = X Aey . 

• We define the differential operators Ax^^ = d"^ and = -|- d'^^. 
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The scaling we consider (we refer to the appendix for physical explanations) leads to 
the study of the scaled Vlasov-Poisson system (for t > 0,x £ := R^/Z^,v G M^): 



dtfe + ^.V,/e + v^^.V^f, + (E, + ^).V,/, = 

= (-v,.^y„-eV,iiK) .... 

-e^A,.,, K - A,.^y, = / Uv - J fdvdx ^ ■ > 

fe,t=0 = fefl > 0, / fefidvdx = 1. 

The quantity f^(t,x,v) is interpreted as the distribution function of the electrons: this 
means that fe{t,x,v)dxdv is the probability of finding particles at time t with position x 
and velocity v; Ve{t,x) and E^(t,x) are respectively the electric potential and force. 

This corresponds to the so-called finite Larmor radius scaling for the Vlasov-Poisson 
equation, which was introduced by Frenod and Sonnendriicker in the mathematical litera- 
ture [9] . The 2D version of the system (obtained when one restricts to the perpendicular 
dynamics) was studied in [9] and more recently in |3] and |10| . A version of the full 3-D 
system describing ions with massless electrons was studied by the author in |17) . In this 
work, we considered that the density of electrons follows a linearized Maxwell-Boltzmann 
law. This means that we studied the following Poisson equation for the electric potential: 

K - e^A^iiF, - A^^V, = J f,dv - J f.dvdx. (1.2) 

In this case it was shown after some filtering that the number density weakly-* converges 
to some solution / to another kinetic system exhibiting the so-called E x B drift in the 
orthogonal plane, but with trivial dynamics in the parallel direction. This last feature 
seems somehow disappointing. 

We observed in |17| that in the case where the Poisson equation reads: 

- e^A^ii^, - A^^V, = J fdv - j fedvdx, (1.3) 

we could expect to make a pressure appear in the limit process e — )• 0, due to the incom- 
pressibility constraint: 

Unfortunately, we were not able to rigorously derive a kinetic limit or even a fiuid 
limit from (ll.ip . This is not only due to technical mathematical difficulties. This is 
related to the existence of instabilities for the Vlasov-Poisson equation, such as the double- 
humped instabilities (see Guo and Strauss |15) ) and their counterpart in the multi- fluid 
Euler equations, such as the two-stream instabilities (see Cordier, Grenier and Guo [8]). 
Such instabilities actually take over in the limit e — )• and the formal limit is false in 
general, unless /^^q does not depend on parallel variables, which corresponds to the 2D 
problem studied by Frenod and Sonnendriicker [9]. 

Actually, we can observe that if on the contrary the initial data f^^ depends only on 
parallel variables, we obtain the one-dimensional quasineutral system: 

dtfe + V\\dxje - dxuVedvJe = 



-e 



^dl^^V, = f f,dv - J f.dvdx^^ (1.4) 

fe,t=0 = fe,0 > 0, / fefldvdx\\ = 1. 

The formal limit is easily obtained, by taking e = 0: 

dtf + v\\ / - 9^11 Vdi,^^ / = 

-e^d^^^V = J fdv - J fdvdxii (1.5) 
ft=o = /o > 0, / fodvdx\\ = 1. 
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In |14| . an explicit example of Grenier shows that the formal limit is false in general, 
because of the double-humped instability: 

Theorem 1.1 ([14]). For any N and s in N, and for any e < 1, there exist for i = 1,2,3,4, 
vfix) e H'{T) with \\vl{x) + 1||h= < e^, \\v^{x) + 1/2\\h- < , \\vl{x) - 1/2||h. < e^, 
||7;|(x) — lilies < e^, such that the solution f^(t,x,v) associated to the initial data defined 
by: 

f^fl{x,v) = 1 for v\{x) <v< V2ix) and vl{x) <v< vl{x) 
= elsewhere. 

We also define /q by: 

fo{x,v) = 1 for -1 <v < -1/2 andl/2 <v <1 
= elsewhere. 

Then f^ does not converge to /q in the following sense: 

liminfsup / \fe{t,x,v) — fo{v)\v'^dvdx > (1-6) 

for any T > and also for T = e", with a < 1/2. 

In order to overcome the effects of these instabilities for the usual quasineutral limit, 
there are two possibilities: 

• One consists in restricting to particular initial profiles chosen in order to be stable 
(this would imply in particular some monotony conditions on the data, such as the 
Penrose condition 1201). 



• The other one consists in considering data with analytic regularity, in which case the 
instabilities ( which are essentially of "Sobolev" nature) do not have any effect. 

Here the situation is worst: by opposition to the usual quasineutral limit (see [6], |14j). 
restricting to stable profiles is not sufficient. This is due to the anisotropy of the problem 
and the dynamics in the perpendicular variables. 

In this paper, we illustrate this phenomenon by formally deriving the following fluid 
system, obtained from the kinetic system (II. ID by considering some physically relevant 
monokinetic data (we refer to the appendix for the detailed formal derivation). 

f dtPe + V±iE^Pe) + 9||(7;||,,/9,) = 

dtV\\^, + V±{E^v\^^J+v\\^,d\\{v\^^J = -ea||(^,(t,x) -5||y,(t,X|,) 

= -V^cf), (1.7) 

-e2(9|^(/>, - Axe/'e = Pe- f Pedx± 

—ed'^^Ve = J Pedx± — 1, 
where: 

• />e(t, x_L, xy) : IR+ X T'^ — 7- can be interpreted as a charge density, 

• ?;||^e(t, x_L, xy) : M"*" x — > M can be interpreted as a "parallel" current density. 

• (peit,x\i) and Ve{t,x) are electric potentials. 
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Although we have considerered monokinetic data, (jl.7p is intrinsicahy a "multi-fluid" 
system, because of the dependence on x±. Hence, we still have to face the two-stream 
instabilities ([8]): because of these, the limit is false in Sobolev regularity and we thus 
decide to study the associated Cauchy problem for analytic data. 

We then prove the limit to a new fluid system which is strictly speaking compressible 
but also somehow "incompressible in average". This rather unusual feature is due to the 
anisotropy of the model. The fluid system is the following (obtained formally by taking 
6 = 0): 

dtp + v^{E^p) + d^\{v^\p) = o 



dtvii + V_l(£^-^?^||) + f||5||(7;||) = -5||p(t, 

= V^Al^ {p- Jpdx^) 
J pdx± = 1. 



Xll 



X8) 



We observe that this system can be interpreted as an infinite system of Euler-type 
equations, coupled together through the "parameter" x±. It has some interesting features: 

• This system is highly anisotropic in x± and xy. The 2D part of the dynamics of the 
equation for p is nothing but the vorticity formulation of 2D incompressible Euler. 
Physically speaking, p should be interpreted here as a density rather than a vorticity. 
The dynamics in the parallel direction is similar to the dynamics of incompressible 
Euler written in velocity. We finally observe that the pressure p only depends on the 
parallel variable and not on x±. 

• This does not strictly speaking describe an incompressible fluid, since (£'-'",U||) is 
not divergence free. Somehow, the fluid is hence compressible. But the constraint 
J pdx± = 1 can be interpreted as a constraint of "incompressibility in average" which 
allows one to recover the pressure law from the other unknowns. Indeed, we easily 
get, thanks to the equation on p: 

c^xii J pv^^dx± = 0. (1.9) 
So by plugging this constraint in the equation on pv^\ : 

dt{pv\i) + V±iE^P\\v\\) + d\\{pvp = -d\\p{t,x\i)p, 
we get the (one-dimensional !) elliptic equation allowing to recover —dx«P- 



9|p(t,X||) = c?|^ / pv'^^dx± 



from which we get: 



5||p(t,X|j) = 9|| j pvfdx±. (1.10) 



• From the point of view of plasma physics, E-^.V± is the so-called electric drift. By 
analogy with the so-called drift-kinetic equations |23) , we can call this system a drift- 
fluid equation. To the best of our knowledge, this is the very first time such a model 
is exhibited in the literature. 

From now on, when there is no risk of confusion, we will sometimes write v and 
instead of v\\ and v\\ , . 
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1.2 Organization of the paper 

The outline of this paper is as fohows. In Section [21 we will state the main results of this 
paper that are: the existence of analytic solutions to (ll.7p locally in time but uniformly 
in e (Theorem I2.ip . the strong convergence to (jl.Sp with a complete description of the 
plasma oscillations (Theorem I2.2p and finally the existence and uniqueness of local analytic 
solutions to (|1.8p . in Proposition 12.11 

Section [3] is devoted to the proof of Theorem 12.11 First we recall some elementary 
features of the analytic spaces we consider (section 13. ip , then we implement an approxima- 
tion scheme for our Cauchy-Kovalesvkaya type existence theorem. The results are based 
on a decomposition of the electric field allowing for a good understanding of the so-called 
plasma waves (section l3.2p . 

In sectionHl we prove Theorem l2.2l by using the uniform in e estimates we have obtained 
in the previous theorem. The proof relies on another decomposition of the electric field, in 
order to exhibit the effects of the plasma waves as e goes to 0. 

Then, in section [5l we discuss the sharpness of our results: 

• In sections 15.11 and 15.21 we discuss the analyticity assumption and explain why we 
can not lower down the regularity to Sobolev. In section [513] we explain why it is not 
possible to obtain global in time results. We obtain these results by considering some 
well-chosen initial data and using results of Brenier on multi-fluid Euler systems [5]. 

• Because of the multi-stream instabilities, studying the limit with the relative entropy 
method is bound to fail. Nevertheless we found it interesting to try to apply the 
method and see at which point things get nasty: this is the object of section 15.41 
where we study a kinetic toy model which retains the main unstable feature of system 

m- 

The two last sections are respectively a short conclusion and an appendix where we 
explain the scaling and the formal derivation of system (II. 7p . 



2 Statement of the results 

In order to prove both the existence of strong solutions to systems (ll.jp and (II. 8p and 
also prove the results of convergence, we follow the construction of Grenier |13| . with some 
modifications adapted to our problem. 

In |13| . Grenier studies the quasineutral limit of the family of coupled Euler- Poisson 
systems: 

dtPe + div(p|,t;^) = 

rot^^ = 
^ ediv£;^ = /^,p^/i(de)-l, 

with (Af, ©,^) a probability space. 

Following the proof of the Cauchy-Kovalevskaya theorem given by Caflisch Grenier 
proved the local existence of analytic functions (with respect to x) uniformly with respect 
to e and then, after filtering the fast oscillations due to the force field, showed the strong 
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convergence to the system: 

dtpe + div(pei'0) = 
dtve + v'q.V {ve) = E 

rot£; = ^ > 

We notice that the class of systems studied by Grenier is close to system (|1.7p . if we 
take X = = x_L and (M, //) = {T'^,dx±), the main difference being that we have to 
deal with a dynamics in Q = x±. 

Hence, we introduce the same spaces of analytic functions as in |13) . but this time 
depending also on Q = x±. 

Definition. Let 6 > 1. We define Bs the space of real functions (f) on such that 

= \mk)\S\'^ < +00, (2.3) 

where T(j){k) is the k-th Fourier coefficient of (p defined by: 

The first theorem proves the existence of local analytic solutions of (ll.7p with a life 
span uniform in e. 

Theorem 2.1. Let 5q > 1. Let Pe(0) and fe(0) he two hounded families of Bsq such that 
J pe{0)dx = 1 and: 



I 



Pe{0)dx± — 1 



< CVe, (2.4) 



then there exists rj > such that for every di g]1,5q[, for any e > 0, there exists a unique 
strong solution {p^,Ve) to / fi. 7\) hounded uniformly in C([0,rj(5o — 5i)[,Bs-^) with initial 
conditions {p^(0),Ve{0)). Moreover, y/ed^^Ve is uniformly hounded in C{[0,rj{6Q — 5i)[, S^^). 

Remark 2.1. • The condition 11 f p^{0)dx± — < C-y/e implies that -y/e9||V^(0) is 

hounded uniformly in B^q (this is the correct scale in view of the energy conservation) . 

• Note that for all t '>0, J p^dx = 1. Hence the Poisson equation —ed'^^V^ = J p^dx± — 1 
can always he solved. 

We can then prove the convergence result: 

Theorem 2.2. Let {p^,Ve) he solutions to the system \1. ?)) for < t < T satisfying for 
some s > 7/2 the following estimate: 

(H) ■■ ^sup (||/>e||H^^,,|| + Mhs^^^^^ + < +00. (2.5) 

Then we get the following convergences 

Pe P, 

i 
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strongly respectively in C{[0,T], H^'^ ^^^) and C{[0,T], H^'^ ]^^^) for all s' < s, and 



strongly in C{[0,T], H^'^^) for all s' < s — 1, and where {p,v) is solution to the asymptotic 
system il.S]) on [0, T] with initial conditions: 

y9(0) = limp,(0), 



v{0) = lini \^v,{0) - J peVedx±{0) 
and X||) are gradient correctors which satisfy the transport equations: 

dtE±+ ( [ pvdxAdx,E± = 0, 



with initial data: 

dx^^E+iO) = lim^a,,, (^-^/^a,||K(0) +i J p,v,dx^iO)^ , 

c>,|l^_(0) = lini^a,.,^ (^-V^c>,.||K(0) -i j p,v,dx^{0)^ . 

As explained in the introduction, due to the two-streams instabihties, we have to restrict 
to data with analytic regularity: the Sobolev version of these results is false in general (see 
[H] and the discussion of Section [S]). 

Remark 2.2. • It is clear that solutions built in Theorem \2.1\ satisfy (H). 

• If instead of (H) we make the stronger assumption, for 5 > 1 

{H') : sup ( llPellfi^ + \\Ve\\Bs + WV^dx V^Wbs) < +00, (2.6) 

t<T,€ ^ ' 

then we get the same strong convergences in C([0, T],^^') for all 5' < 6. 

Using Lemma l3.il {ii),{iv), the proof under assumption {H') is the same as under 
assumption (H). 

• The "well-prepared" case corresponds to the case when: 



limdx,, / PeV^dx±{0) = 0. 
e^O II J 



Then there is no corrector. 



With the same method used for Theorem 12.11 we can also prove a theorem of existence 
and uniqueness of analytic solutions to system (jl.Sp . 
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Proposition 2.1. Let 5o > 6i > 1. For initial data p(0),f(0) G Bs^ satisfying 

p(0) > 0, (2.7) 

p{0)dx± = 1 (2.8) 



and 

9|| y p(0)w(0)dx_L = 0, (2.9) 

there exists rj > depending on 5o and on the initial conditions only such that there is a 
unique strong solution (p,U||,p) to the system il.8\) with p,v £ C([0,ry(5o — 6i)[,Bs) for all 
6 < 6i. 

3 Proof of Theorem EH 

3.1 Functional analysis on Bs spaces 

First we define the time dependent analytic spaces we will work with. 

Let /3 be an arbitrary constant in ]0, 1[ (take for instance /? = 1/2 to fix ideas) and 
1] > a parameter to be chosen later. 

Definition. Let 6o > I. We define the space B^^ = {u G C^{[0,r]{6o — 1)], Bsg^t/ri)}) 
endowed with the norm 

\\u\\s,= sup (\u{t)\s + {So-S--)^\Vu{t)\s) 

l<S<5o ^ ^ 

0<t< r/((5o - S) 



where the norm \u\s was defined in l[2.3]) : 

\u\s = ^ \Tu{k)\5^''^ < +00, 

We now gather from |13) a few elementary properties of these spaces, that we recall for 
the reader's convenience. 

Lemma 3.1. For all S > 1 

(i) The spaces Bs and B^ are Banach algebra. 

(a) If 5' < 5 then Bg C Bs' , the embedding being continuous and compact. 
(Hi) For all s G M, i?^ C H'^ , the embedding being continuous and compact, 
(iv) For all 1 < 5' < 5, if (p e Bs, 

(v) If u is in B^^^ and if 5 + t/rj < Sq then 



\dl^^u{t)\s<2'+^\\u\\sM^o-S-l^)-^-\ 
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For further properties of these spaces we refer to the recent work of Mouhot and Villani 
|19j . in which similar analytic spaces (and more sophisticated versions) are considered. 

Proof. We give an elementary proof for [ii) which is not given in |13) . The embedding is 
obvious. We consider for iV G N the map ztv defined by: 

\k\<N 

We then compute: 

\{id-tN)<pw = E \^m\s'^'^ <[j) E \mk)\s^''^ <[j) \4>\s. 

\k\>N ^ ' \k\>N ^ ' 

So the embedding Bf, C -B^/ is compact as the limit of finite rank operators. 

- 5- ■ 
2 

We will also need the following elementary observation: 
Remark 3.1. Let (j) ^ B^. Then: 



For (-y), take 5' = 5 + 1 and apply {iv). We refer to |13| for the other proofs. □ 



(j)dx± 



S 



Proof. We simply compute: 



□ 



3.2 Description of plasma oscillations 

To simplify notations, we set i?^ y = —dx^^Ve{t,x\\) (which has nothing to do with E^). In 
this paragraph, we want to understand the oscillatory behaviour of We will see that 
the dynamics in x± does not interfer too much with the equations on -E^ji, so that we get 
almost the same description of oscillations as in Grenier's paper |13) . 

First we differentiate twice with respect to time the Poisson equation satisfied by V^: 

ed^dx^^E,^\\ = J p,dx±. (3.1) 
We use the equation on to compute the right hand side of (13. ID . 

dt j pedx± = - j \7±{E^pe)dx±-dx^^ j p^v^dx^. (3.2) 

^ V ' 

=0 

Then we integrate with respect to x± the equation satisfied by p^v^, that is: 

dtiPeVe) + y±{E^PeVe) +dx^^{v^Pe) = - Pe{edx^^(t)e{t, x) + da:^^V^{t, X\\)) 

and we get: 

- dt J PeVedx± = j Pev'^dx± + || J p^dx± - j p^{ed:c^^4>e)dx±, (3.3) 
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so that: 



Thus it comes: 



n '-'X 



j Pevldx^ + ed:^^^[E^^\\dx^^E^\\ - d^^^ j peiedr,^^(j)^)dx±. (3.4) 



Equation (|3.4p is the wave equation allowing to describe the essential oscillations. At 
least formally, this equation indicates that there are time oscillations with frequency 

and magnitude created by the right-hand side of the equation which acts like a source. 
We observe here that the source is expected to be of order 0(1): indeed, by assumption 
on the data at t = 0, we can check that this quantity is bounded in a space. 

In particular if we want to prove strong convergence results we will have to introduce 
non-trivial correctors in order to get rid of these oscillations. We notice also that p.4p is 
very similar to the wave equation obtained in |13| (the only difference is a new term in the 
source), so that most of the calculations and estimates on E^ y we will need are done in 

m 



3.3 A priori estimates 

We have just ob 
average in time 



We have just observed that iJ^ y roughly behaves like -^e^**/^. Hence if we consider the 

ft 

a= / ^,,|i(s,X||)ds, (3.5) 
Jo 

we expect that is bounded uniformly with respect to e in some functional space. We 
also introduce the translated current (which corresponds to some filtering of the time 
oscillations created by the electric field): 

We=Ve-Ge, (3.6) 

so that system (|1.7p now writes: 

dtPe + '^±{E^Pe) + 5||((u;e + G,)p,) =0 

dtw, + V±{E^{w, + Ge)) + {w, + Gja,, {v, + a)) = -ed\\Ut, x,, ). ^'"''^ 

The goal is now to prove some a priori estimates for G^, Pe and w^- We are also able to 
get similar estimates on Ej- and edxt^^<pei thanks to the Poisson equation satisfied by (f)^. 

3.3.1 Estimate on G^ and \/e£^e,|| 

We use Duhamel's formula for the wave equation (|3.4p to get the following identity: 

rt 



F\\G,{t,k\\) = l_cos(^) T\\g,{s,k\\)^ds + T\\G'i, 



(3.5 



denoting by the Fourier transform with respect to the parallel variable only and the 
Fourier variable and where: 

= <9^|| j Pevldx^ + e5^|| - 5^11 j Pe{edx^^(l)e)dx±, 
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= (0, X|| ) sin (^-^^ - edtE.^^i (0, X|j ) (^cos (^-^^ - 1^ . 

We now estimate HGeH^Q. 



(3.9) 



\Ge\5 < 



1 



,t- s 



II \ik 



[1 - cos(— ^)]J"||5e(s,A;|i) 



+ |G°| 



Thanks to Remark 13. II and Lemma l3.ll . (i): 
Similarly, we prove: 



< 



(3.10) 



<kSo-6-'-)-P\\V~eE,,\\\\l, 
5 2 rj " " 



j dx^^ {p^{edx^^4)e)^ dx± 



Thus, we have: 

ft 



\G.\s <C j{5,-6- -)^~^\\\pMW\\1 + \\VeE,,\\\\l + MsJedx.Us,) + \G%. 







Likewise, one can show (this time we use lemma l3?T| {v)): 

|9x„a|5 < C fid, - 5 - -Y~^-'H\\Pe\\Ml + \\V~eE,,\\\\i + WPehohd^Aho) + \d.,G%. 
Jo V 

Hence using the elementary estimates 

ds 



I 

Jo 



^ 2 ;:l-/3 
< ■>]-- 



f 

Jo 



{So-6-^r - '1-/3 



ds 



'0 ' 



^ 2r/ t 



and \\ve\\sa < WweWso + ||G,||5(,, we get: 



liaiUo < vCiSo,P) [iMo + WGeWsofMso + WVeE^ll + \\pe\\5o\\^d,,ct>e\\5o) + \\G%o. 

(3.11) 

If we compare two solutions {w^^\ p^^^) and {w^'^\ p^'^^) with the same inital data we obtain: 



xiWwPWs, + H'^Ws, + ||G«IUo + l|G?)|Uo)(llp^')|Uo + \\P^'%o) 
+ (ll^^^lll. + H'X + + ||Gf)|||,)(||p« 

+ iiv^i?;^ - vi^;^?ii,„(iivi^!^)|i,„ + ii^/^i?;')|Uo)). (3.12) 
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Likewise we get the same estimates on ||y^£'^ y since we have the formula: 

= (^_L[sin(^)]J^,5,(5,A:i|)) + (3.13) 

with 

El. =E,|,(0,x)cos(^) + ^/iatS,|j(0,x)sin(^). (3.14) 
3.3.2 Estimate on Ej- and edx^^4>t 

We now use the scaled Poisson equation satisfied by (/j^ to get some a priori estimates. 

The principle here is to look at the symbols of the operators involved in the Poisson 
equations. Accordingly, we compute in Fourier variables: 

e^k^Fcl)^ + |fc_L|^J"(/>, = J=-(^p,-j pedxj_^ . (3.15) 

Thus it comes: 

Since J(p^ — J p^dx±)dx± = 0, we have for all /cy: 



J" 

Thus it comes, for all A; i , ki\ : 



(^Pe - J Pedx±^ (0, A:||) = 0. 



\J^{Pe - f Pedx±)\ 



|-F<Ae| < ' ^^'J, (3.16) 

In particular we easily get, using the relation E^ = —V'^cpe- 

\J^E-\ < \^^P^-f^P^d^-^\ < 1^ -J.^dx, 
Hence: 

\\E^\\5o<CMso. (3.17) 
Likewise, since ab < i(a2 + ^2) and |A;_l| > 1: 

and consequently: 

\\ed,^^<l)e\Uo<C\\Pe\\5o- (3.18) 
And if we compare two solutions with the same initial data: 

\\E^'^'^ - + Ik9.„# - ^5.„#||,, < - pi'^so- (3.19) 
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3.3.3 Estimate on and 

We now use the conservation laws satisfied by and to get the appropriate estimates. 
The density p^ satisfies the equation: 

dtp, + V^{Etpe) + + = 0. 

With the same kind of computations as before and thanks to estimate (|3.17p we get: 



PeU < / \dtPe\s + \PeiO)\s 

Jo 



/■* S 

< l|Pe(0)ll<5o + / ('^'O-'^ )~^ll/Oj5o(ll/5e||5o + IkelUo + ll'^elUo)'^^- 

Jo V 



Similarly we estimate {dx^Pels by differentiating with respect to Xi the equation satisfied 
by p^. Finally we get: 

WpeWso < vCWpehoiWpeho + lkJ<5o + llCelUo)- (3-20) 

If we compare two solutions with the same initial conditions, we get likewise: 

X iW^-p^^^Wso)). (3.21) 
In the same fashion, we estimate the 6o norm of w^: 

\\w,\\so < r]C {{\\w,\\5, + l)\\Pe\\so + ilMso + WGeWso? + I|e5|,0e||5o) > (3-22) 
and if we compare two solutions with the same initial data: 

+ i\\4'^-^i'^\\So + \\pP-pP\\s,) 

+ iH'^\\So + H'^\\So + \\Gi'^\\So + \\G^^^\\So) 

||e9||0W-e5||#||5o). (3.23) 
3.4 Approximation scheme 

We use the usual approximation scheme for Cauchy-Kovalevskaya type of results (|7]). 
We define w'^, G", V^", (p^ by recursion: 
Initialization For < t < ri(6o — 1), we define: 

p'At)=Pe{0), 

w^, = v,{0) - GO, 
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J p^dx. Finally, is given by formula G° = - £ 9|| V^^. 



and -e^a^^y^o = pO- 
Recursion We define /o""^^, zu""^^ by the equations: 

f dtp^+^ + V^iE^'^'p^) + 9,1 (« + = 

1 9t<+i + Vx(^e^'"« + G^)) + « + G^)dii « + G^)) 



-ea|0^(t,2;|| 



with the initial conditions: /5"+^(0) = Pe{0) and 11;"+^ = ^^(0) — G°. 
Then we can define 0"+^ as the solution to the Poisson equation: 



2fl2 in+1 



A, 



.,n+l 



-V 



Similarly, 



,n+l 



dx. 



Then we can define G""*"^ with the formula: G"+^ = — Jg Syy^""'"^. 

Now let Gi be a constant larger than ||/3e(0)||5p, ||w;e(0)||5(,, ||Ge(0)||5,j, ||-v/e£'e(0)||5Q and 
all the other constants in the previous estimates. It is possible to choose rj small enough 
with respect to Gi to propagate the following estimates by recursion (we refer to [13] for 
details; we use in particular estimates p.l2p . p.2ip . ()3.23p ). There exists G2 > Gi, for all 
n > 1: 



(i) 



WPeWSo < G2, 
ll<ll5o < C2, 
||G"||5g < G2, 

\\VeE%\\5o<C2. 



(ii) 



\\w- 
\\G"' 



n-l|| ^ C2 
re 1 1 00 — 2" ' 
n-l|| < C2 
^e 1 1 00 — 2" ' 
nn-l\\ < C2_ 
1 1 00 — 2" ' 



^,11 



n.-li 



l<5o 



< 



C2 



This proves that the sequences p'^,w^,G'^,^E^,E^ '^^edx^^cj)^ are Cauchy sequences 
(with respect to n) in B^^, and consequently converge strongly in i?^^, the estimates being 
uniform in e. It is clear that the limit satisfies System (|1.7p . 

The requirement 61 < 6q and the explicit life span in Theorem 12.11 come directly from 
the definition of the B^^ spaces. 

For the uniqueness part, one can simply notice that the estimates we have shown allow 
us to prove that the application F defined by: 



E{Pe,We) 



jt{-V^{Etpe) - d\\{{w, + G,)p,)))ds ' 
jl{-V^{E^{w, + G,)) - [w, + G,)d\\{v, + G,)) - e9||0e(t,x||))ds 



is a contraction on the closed subset B of x Bs^ , defined by: 

B = {p,w e Bso', WpWso < G, \\w\\so < G} , 

with G large enough, provided that r] is chosen small enough. The uniqueness of the 
analytic solution then follows. 
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Proof of Proposition 12.11 

We can lead the same analysis as for the proof of Theorem 12. H but even simpler since here 
we do not have to deal anymore with the fast oscillations in time. The only slightly different 
point is to estimate the norm of —d\\pds = J pv^dx±ds, which is straightforward: 

t 

d\\pds < vCWpWsoMIo- 

So 

Then as before, we can use a contraction argument to prove the proposition. 



4 Proof of Theorem 12.2 



Step 1: Another average in time for 

We have observed previously that the wave equation (|3.4p describing the time oscillations 
of II was the same as the one appearing in Grenier's work, except for a slight change in 
the source. Therefore the following decomposition taken from |13| identically holds: 

Lemma 4.1. Under assumption (H), there exist vector fields El,E'^ and such that 
E'^ II = El + E1 and a positive constant C independent of e such as: 

(i) II ^^,^11^^(^.-1) < C. 

(ii) dtWe = El, ||W,||^^(^.-i) <C andW^^O in L^. 

(ttt) W'iO) = -edtE.^iiiO) = J p,iO)v,{0)dx^. 

(iv) 11^2 11^^ < c. 

(v) J Eldx\\ = J Eldx\\ = 0. 

Idea of the proof. The idea in order to build E^ is to cut off the essential temporal oscilla- 
tions (of frequency -^). Hence, we can define E^ defined by its Fourier transform: 

Jll^2(i^^lj)^__ ^11^,^11 (s,A;|j)ds 

and El = y — E'^, so that El gathers the essential information on the oscillations. 
Let us refer to |13| for details. 

□ 

Step 2: Uniform bound on E-^ and 9^ 

Under hypothesis {H) we clearly get that E^ and dx^^(t>e are bounded in L'l°{H^~^) uni- 
formly with respect to e (we do not need any gain of elliptic regularity). 
Since 

JiP^~ j Pedx±)dx± = 0, 

we easily check that: 

||f/'e|k5^,,|l ^\\P- j P^^^±IIh|^..||- 

Hence the result. 
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Step 3: Passage to the strong limit 

Let We = Ve — We- According to Lemma |4.H We is uniformly bomided in L'^{H^~^). On 
the other hand, we have : 

dtw, + V±{E^We) + wA^^w, = -e9,||</)e + - w,d,^^W, - W,d,^^w, - W,d,^^W,. (4.1) 

(Notice that V±{E^We) = WeV±{E^) = 0.) 

Thus, using the uniform bounds, we can see that dtw^ is bounded in L^(i?*~^) and 
thanks to the Aubin-Lions lemma (see J. Simon [21]), converges strongly (up to a 
subsequence) to some function w in C{[0,T],H''-^) for all s' < s. 

According to Step 2, edx^^4'e ^ in the distributional sense. The following convergence 
also holds in the sense of distributions, according to Lemma l4.ll 

and Wedx^^W^ + E'^ weakly converges to some function F since it is clearly bounded in 
Furthermore, since: 

J (w,dx^^W, + E^'^ =0, 

this implies that J Fdx^^ = 0, and thus there exists p such that F = —dx^^P- 

Since E^ is uniformly bounded in L'^{H'^~^), it also weakly-* converges, up to a 

subsequence, to some function E-^. 

We now use the strong limit of We in C{[0,T], H^'^^) in order to pass to the limit in 

the sense of distributions in the convection terms. As a consequence, it comes, passing to 

the limit in the sense of distributions: 

dtw + \/±{E^w) + wdx^^w = -dx^^p. (4.2) 
The equation satisfied by is: 

The proof is similar for which converges strongly, up to a subsequence, to p in 
C{[0,T],H' ) for aU s' < s. One can likewise take limits in the Poisson equations. We 
finally obtain (II. 8p . By uniqueness of the solutions to (II. 8p , the limits actually hold without 
extraction. 

Step 4: Equations for the correctors 

The final step relies on the following lemma proved in Grenier's paper |13) (the main point 
is to notice that the application / i— t- e^**/^/ is an isometry on L°°{H^) for any s.) 

Lemma 4.2. There exist two correctors i?+(t,a;|j) and i?„(t,X||) in C{H^~^) such that, for 
all s' < s: 

• \\V~eEl - e^*/v^S+ - e-^*/v^S_||c(H.'-i) ^ 0, 
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In particular we can deduce that: 



(and similarly e'■*/'^^/eEl E_). 

Then, the idea is to use Lemmas 14.11 and 14.21 and the wave equation (j3.4p in order to 
obtain the equations satisfied by E± . By elementary (but rather tedious) computations we 
get: 

9t(5,||^±) + (^j pvdxj}j a,||(a,||^±) = 0. 

Lemma [4. 21 also provides the initial conditions for E±. 
The proof of the theorem is now complete. 



5 Discussion on the sharpness of the results 
5.1 On the analytic regularity 

Let us recall that the multi-fiuid system (12. 2p is ill-posed in Sobolev spaces, because of 
the two-stream instabilities (remind that this is due to the coupling between the different 
phases of the fluid) . 

For system (II. Sh . we expect the situation to be similar. Due to the dependence on x± 
and the constraint J pdx± = 1, system (jl.Sp is by nature a multi-fluid system. Neverthless, 
one could maybe imagine that the dynamics in the x± variable could yield some mixing in 
x± and X|| (in the spirit of hypoellipticity results) and thus could perhaps bring stability. 
Here we explain why this is not the case. 

The idea is to consider for (II. 8p shear flows initial data. This will allow to exactly 
recover the multi-fluid equations (|2.2p . We take: 

E^ = (0,v9(xi,rE||),0), 

and consequently po = V± A Eq = —ip'{xi,u). We also assume that vo{xi,x\^) does not 
depend on X2- 

Then we observe that: 

^±{E^Po) = 0, 
V^iE^vo) = 0. 

With such initial data, system (jl.Sp reduces to: 

5tP + 9||(u||p) = 

dtv\\ +v\\d\\{v\i) = -d\\p{t,x\\) (5.1) 
/ pdxi = 1, 

and we observe that there is no more dynamics in the x± variable. This is nothing but 
system (12. 2p in dimension 1, with M = [0, 1[ and fi the Lebesgue measure. 

Now, let us consider measure type of data in the xi variable for p and v (this corresponds 
to a "degenerate" version of the shear flows deflned above). In particular if we choose: 

1 / ^ 1 



17 



we get: 



\ \ (5.2) 

^0 = 2^xi=i'^o,i(a;|j) + 2'^xi=i^o,2(x|j) 

and we obtain the following system for a = 1,2: 

QtPa + Q\\ [VaPa) = 

dtVa + Vad\\ [Vo,) = -d\\p{t, X|j ) (5.3) 

Pl + P1 = 1- 

This particular system was given as an example by Brenier in [5] to illustrate ill-posedness 
in Sobolev spaces of the multi-fluid equations. 

We denote q = piVi. Using the constraint + = 1, we easily obtain that 

2 n 1 
P\\=-1 — + 

We can then observe that the system: 

f dtpi + diiq = 



5t9 + 9||(f^) = -pi9||p(t, 



(5.4) 



is elliptic in space-time, and consequently it is ill-posed in Sobolev spaces. 

Actually this example is not completely satisfying, since it is singular in xi. Nevertheless 
we can consider the convolution of this initial data with a standard mollifier, which yields 
the same qualitative behaviour. 



5.2 On the analytic regularity in the perpendicular variable 

We observe that if the initial datum (p(0),u(0)) does not depend on xy, then the fluid 
system (II. 8p reduces to: 

f dtp + V^{E^p) = 
^ J pdx± = 1. 

Thus, p satisfies 2D incompressible Euler system, written in vorticity formulation. This 
systems admits a unique global strong solution provided that p(0) G H^(T'^) (with s > 1), 
by a classical result of Kato |18| and even a unique global weak solution provided that 
p(0) G L°°(T2), by a classical result of Yudovic [21]. 

In the other hand, satisfied a transport equation with the force field E-^. If we 
assume for instance that vq is a bounded measure, then using the classical log-Lipschitz 
estimate on E^, we get a unique global weak solution v^i^ by the method of characteristics. 

One could think that it should be possible to build solutions to the final fiuid system 
(II. 8p with similar "weak" regularity in the x± variable (while keeping analyticity in the 
X|| variable). Actually this is not possible in general: this is related to the fact that E-^ 
depends also on xy and this entails that we also need analytic regularity in the x± variable 
to get analytic regularity in the xn variable (see estimations such as p.20p ). 
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5.3 On the local in time existence 



In [5], Brenier considers potential velocity fields, that is velocity fields of the form v@ = 
Vx^e, for the multi-fiuid system: 

f G = 1,...,M M G N* 

dtpe + div(/9e-ue) = 

dtve + ve.V{ve) = E (5.6) 
rot £■ = 

I 1.0=1 Pe = 1- 
In this case the equation on the velocities becomes: 

5t$e + ^|V,$e|2+p = 0. (5.7) 

It is proved in [5] that any strong solution satisfying 

inf pe{t, x) > 

0,t,x 

can not be global in time unless the initial energy vanishes: 

M 



^ J pe,t=ohe,t=opd2; = 0. (5.8) 



e=i 

This striking result relies on a variational interpretation of these Euler equations. Using 
the same particular initial data as in section 15.1) this indicates that for system (II. 8p also, 
there is no global strong solution, unless there is no dependence on x± or xy. 

We observe that if the initial datum (p(0),u(0)) does not depend on x±, the fluid 
system (jl.Sp does not make sense anymore (as for incompressible Euler in dimension 1). 
When the initial datum (p{0),v{0)) does not depend on xy, we have seen that we recover 
2D incompressible Euler and there is indeed global existence (of strong or weak solutions). 

5.4 The relative entropy method applied to a toy model : failure of the 
multi-current limit 

It seems very appealing to try to use the relative entropy method (which was introduced 
by Brenier [1] for Vlasov type of systems) to study the limit, as it would open the way to 
the study of the limit for solutions to the initial system (II. ID with low regularity. The only 
requirement would be that the two first moments of the initial data for (|1.1|) are in a small 
neighborhood (say in topology) of the smooth initial data for the limit system (II. Sp . 
Nevertheless it is not possible to overcome the two-stream instabilities in this framework. 
We intend to show why. 

Let us consider the toy model: 

dtf! + v.V.f^ + E,.V,f! = 
E, = -V.K 

-eA,K = ///>d/i-l ^ ' 



m = 0) = /, 



e,0- 



with t > 0, X G T'^, u G M'^ and where 9 lies in [0, 1] equipped with a positive measure ji 
which is: 
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either a sum of Dirac masses with total mass 1, such as: 

7V-1 



1=0 

In this case, we model a plasma made of N phases. 

• or the Lebesgue measure, in which case we model a continuum of phases. 

Actually, we could have considered more general Borel measures but we restrict to these 
cases for simplicity. This system can be seen as the kinetic counterpart of a simplified 
version of (II. 7p . which focuses on the unstable feature of the system. Of course we could 
have considered directly the fluid version, that is: 

dtpl + vMO = 

df.+<.Vul = E, (5.10) 
-eA^K = / / f.dvdp - 1 

but the proofs are essentially the same and the study of system (j5.9p has some interests of 
its own. 

We consider global weak solutions to (j5.9p , in the sense of Arsenev [1] . We recall that 
the energy associated to ()5.9p is the following non-increasing functional: 



£e{t) = \j j fM''dvdxdp+]^e j |V,.K|'dx. (5.11) 

We assume that there exists a constant K > Q independent of e, such as <?e(0) < K. 
This implies that for any e and t > 0: 

8,{t) < K. (5.12) 
Let {p^,u^) be the local strong solution, to the system: 

V + v..(pV) = 

dtu^ + u^.V^u'^ = -V.^V (5.13) 
/p% = l. 

with inital data (pq,Uq) (which we a priori have to take with analytic regularity). The 
"incompressibility" constraint reads: 



. J p\^dp = 0. (5.14) 



Vx- j p\^dp = 0. 

Following the approach of Brenier [3] for the quasineutral limit with a single phase, we 
consider the relative entropy (built as a modulation of the energy 

n,{t) = \jj fe\v - u^{t,x)\^dvdxdp + ^ej |V,K - V,V\^dx. (5.15) 

We assume that the system is well prepared in the sense that ^^^(0) — t- 0. The goal is 
to find some stability inequality in order to show that we also have 7ie{t) — t- 0. 

Our aim here, is to show why the method fails unless actually does not depend on 
0. This can be interpreted as the effect of the two-stream instabilities [8]. 
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We have, since the energy is non-increasing: 



^^.(i) < j j dtf! QhT - y-u^^ dvdxdii + j j f^dt QkT - v.u'^^ dvdxdfi 

J dt\V,V\^dx-e J V,V,.dtV^Vdx - e J dtV,V,.V,Vdx. 

(5.16) 

We clearly have e J dtlVxVl'^dx = 0{e). Moreover, we get, using the conservation of 
energy. 



/ 



Vr^Ve.dt^xVdx 



<V^\\V^v,v,\\L2\\dtV,v\\L2, 



which is of order 0{^/e). 

For the last term of (I5.16p . we compute: 



-e / dtV^Ve.VxVdx =e j dt^^V.Vdx 

= - e I dt{ I f^dvdn ] Vdx 
= + 



y" V,. (^y" f^vdvdi?j Vdx 



(5.17) 



f^vdvdfij V.j:Vdx. 
In this computation we have used the local conservation of mass: 

dt j f^dv + V^. {^j vf!dv^ =0. 
In the other hand we can compute: 



f^{u^ -v).{u^ -v).Va:Uedvdxdfi + j J f^{u^ - v).{dtu^ + u^Va:U^)dvdxdn 
f^E^.u^dvdxdfi. 

(5.18) 

All the trouble comes from this last term. When no assumption is made on , it can be 
of order 0{l/y/e). This wild term can be interpreted as the appearance of the two-stream 
instabilities. 

Therefore we have to make an additional assumption in order to avoid this instability. 
This is done by assuming that initially does not depend on 9 (which yields that does 
not depend on by uniqueness , in which case we can write: 

ue = u 

and consequently, we have 

- J f^E^.udvdxdfi = J (eA^V; - 1) E^.udx. (5.19) 
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In addition, the incompressibihty constraint becomes Vx-U = 0, and thus: 



EfUdx = / VfV x-udx = 0. 



Furthermore, we have: 



j (^j f^dvdp^ u.VxVdv = j u.VxV-e j /^^VeU-V (5.20) 

The first term is equal to according to the incompressibihty constraint, while the second 
is of order O(y^), by the energy inequality. 
We finally get the stability inequality: 



+ 



Ueit) <nM + RM + c [ \\Vxu\\ne{s)ds 

Jo 



(5.21) 



with Re{t) — )• as e goes to and the last term is by definition of n, V. 

As as result, by Gronwall's inequality, T-Le{t) — )• 0, uniformly locally in time. 

To conclude, we notice that := f ffdv is uniformly bounded in L'^{L\^). By the 
energy inequality, we can easily show that this is also the case for Jf := / ffvdv. Thus, 
up to a subsequence, there exist pg and such that weakly-* converges in the sense of 
measures to p^ (resp. to J^). Passing to the limit in the local conservation of charge, 
which reads: 

dtpl + V,. jf = 0, 

we obtain: 

dtP^ + V,/ = 0. 

The goal is now to prove that = p^u. 

By a simple use of Cauchy-Schwarz inequality, we have: 



// 



. "^'^ dxdp < I I f^\v-ufdvdxdp. (5.22) 



Using a classical convexity argument due to Brenier [6], the functional {p, J) i— )■ J ^^"""^^ dxdp 
is lower semi-continuous with respect to the weak convergence of measures. We finally ob- 
tain by passing to the limit that: 

,f = p'u. 

By uniqueness of the solution to the limit system, provided that the whole sequence p^ g 
weakly converges to Pq, we obtain the convergences without having to extract subsequences. 
Finally we have proved the result: 

Proposition 5.1. Let {f^,E^) be a global weak solution to l[5.y\) such that the local conser- 
vation of charge and current are satisfied. Assume that for some smooth functions (pq,uo) 
(we emphasize on the fact that uq does not depend on 0, in order to avoid two-stream 
instabilities ) satisfying: 

i^'^^ = 1' (5.23) 
Vx-uo = 0, 
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we have: 



\ j j fe{t = Of\v - uo{x)\'^dvdxdn + j WxVe{t = 0) - V^V\^dx (5.24) 
and J ffdv pg in the weak sense of measures. Then, 

\jj f!\v-u{t,x)\'^dvdxdp+]^e j |V^K- V^y|2(ix^0, (5.25) 



2 

where {u, V) is the local strong solution to the incompressible Euler system: 

dtu + u.VxU = —VxV 

S/rU = 0. 



(5.26) 



Moreover, := J f^dv converges in the weak sense of measures to p^ solution to: 

dtp^ + n.V^./ = 0. (5.27) 
and := J f^vdv converges in the weak sense of measures to p^u. 

6 Conclusion 

In this work, we have provided a first analysis of the mathematical properties of the three- 
dimensional finite Larmor radius approximation (FLR) , for electrons in a fixed background 
of ions. We have shown that the limit is unstable in the sense that we have to restrict to 
data with both particular profiles and analytic data. In particular we have pointed out that 
the analytic assumption is not only a mathematical technical assumption, but is necessary 
to have strong solutions. In addition the results are only local-in-time. 

On the other hand, we proved in |17| that the FLR approximation for ions with massless 
electrons is by opposition very stable, in the sense that we can deal with initial data with 
no prescribed profile and weak (that is in a Lebesgue space) regularity. 

This rigorously justifies why physicists rather consider the equations on ions rather 
than those on electrons, especially for numerical experiments (we refer for instance to 
Grandgirard et al. [T2]). 



7 Appendix : Formal derivation of the drift-fluid problem 
Scaling of the Vlasov equation 

Let us recall that our purpose is to describe the behaviour of a gas of electrons in a neu- 
tralizing background of ions at thermodynamic equilibrium, submitted to a large magnetic 
field. For simplicity, we consider a magnetic field with a fixed direction ey (also denoted 
by Cz) and a fixed large magnitude B. 

Because of the strong magnetic field, the dynamics of particles in the parallel direction 
ey is completely different to their dynamics in the orthogonal plane. We therfore consider 
this time anisotropic characteristic spatial lengths: 

X± Xll 

x± = — , x± = — 
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f{t,x±,xpv) = ff{i,x±,x\^,v) V{t,x±,x^^) = VV{i,x±,x\\) E{t,x±,xi\) = EE{i,x±,x\\). 

The subscript _L stands for the orthogonal projection on the perpendicular plane (to 
the magnetic field), while the subscript || stands for the projection on the parallel direction. 
This yields: 



P^ = (-itV,,v;,-^v,„k 
-ff A.^v; - f[A,„v; = efvl [Jhdv - 1 

. fe,\i=0 = f0,e, fL^vf,^ J fo,edvdx = 1. 



(7.1) 



We set = ^ : this is the cyclotron frequency (also referred to as the gyrofrequency) 
We also consider the so-called electron Larmor radius (or electron gyroradius) defined 
by: 

Vth rnvth „s 

This quantity can be physically understood as the typical radius of the helix around axis 
ey described by the particles, due to the intense magnetic field. 
The Vlasov equation now reads: 

dj, + ^nrv^.Vs^J, + ^J]T^||.Vi|| A + (-g^f^r^e + ^tv A e,,^ .Vj, = 0. 

The gyrokinetic ordering consists in: 

ilT = -, ^ = e. 

The spatial scaling we perform is the so-called finite Larmor radius scaling (see Frenod 
and Sonnendrucker [9] for a reference in the mathematical literature) : basically the idea is 
to consider the typical perpendicular spatial length L± with the same order as the so-called 
electron Larmor radius. 

On the contrary, the parallel observation length is taken much larger: 

This is typically an anisotropic situation. 

This particular scaling allows, at least in a formal sense, to observe more precise effects 
in the orthogonal plane than with the isotropic scaling (studied for instance in |11|): 

In particular we wish to observe the so-called electric drift E-^ (also referred to as the 
E X B drift) whose effect is of great concern in tokamak physics (see [16] for instance). 
The quasineutral ordering we adopt is the following: 

^ = V-e. (7.4) 
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After straightforward calculations (we refer to [9] for details), we get the following 
Vlasov-Poisson system in dimensionless form, for t > 0,x = (x_l,X||) G x T, u = 
{v±,vu) G ]R2 X M: 







^e = (-V,^K,-eV.||K) 
-e^A^llK - A,^V, = J f,dv - J hdvdx 

. /e,t=0 = fefi- 



Remark 7.1. It seems physically relevant to consider scalings such as: 



(7.5) 



(7.6) 



with a > 1. However with such a scaling, the systems seem too degenerate with respect to 
€ and we have not been able to handle this situation. The scaling we study is nevertheless 
relevant for some extreme magnetic regimes in tokamaks . 



Hydrodynamic equations 

In order to isolate this quasineutral problem, thanks to the linearity of the Poisson equation, 
we separate the electric field in two parts: 

' E = El + El 

El = (-V,.^V;,-eV,||K), 
< -e2A,||V;i - A^^y/ = / f,dv - f hdvdx^, (7.7) 

—e^xitye = / fedvdx± — f fedvdx. 

In order to make the fast oscillations in time due to the term ^.V^- disappear, we 
perform the same change of variables as in |10) . to get the so-called gyro-coordinates: 

Xg = X± + V^,Vg = v±. (7.8) 

We easily compute the equation satisfied by the new distribution function ^^(t, Xg,Vg, uy) = 
fe{t,x,v). 

dtQe + v\\dx^^ge + El^^^{t,Xg - Vg)dv^^ge + E'^{t,Xg^ii)dy^^g^ 
+El±{t,xg - vf).{V,^g, - V^^J + ^^'^V„^g, = 0. 

Notice here that in the process, the so-called electric drift E-^ appears since: 
-El^it,xg - vj).V^^g, = El'^{t,xg - vf).V,^g,. 
The equation satisfied by the charge density Pe = J gedv states: 

dtPe + / v\\gdv + V^^ / El^{t, Xg - vf)g,dv = 0, (7.9) 
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and the one satisfied by the current density Je = J Qevdv ^= ^ J^*^^"*"^^^^ following: 



We now assume that we deal with special monokinetic data of the form: 

9e{t,X,v) = p^(t,2;)lj,||=„||^^(t^^)l„^=0- (7.11) 

This assumption is nothing but the classical "cold plasma" approximation together with 
the assumption that the transverse particle velocities are isotropically distributed (which 
is physically relevant, see |22) ) : in other words, the average motion of particles in the 
perpendicular plane is only due to the advection by the electric drift E-^. 

For the sake of readability, we denote by now V^g = V_l and V^y = Vy. Then we get 
formally the hydrodynamic model: 



f dtPe + V±{E^Pe) + d\\ {V\\^,p,) = 

dt{PeV\le) + ^ ±iE^PeVle) + dl\{peV'^^J = -e9|| (t, x)pe " <9|j K (t, )Pe 

E^ = -V^cj), 

-e^dj(l)e - /^^4>e = Pe- f Pedx± 

-edf^Ve = Jpedx± - 1 



(7.12) 



One can use the first equation to simplify the second one (the systems are equivalent 
provided that we work with regular solutions and that p^ > 0): 



f dtp, + V±{E^Pe) + d\\{v\\^,p,) = 

dtv\\^, + V±{E^v\\^,) + v\\^,d\\{v\i^,) = -e5||0,(t,x) -5||ye(i,a;||) 
E^ = -V^6, 



A 



Pe- J Pedxj 



(7.13) 



-edf^Ve = J Pedx± - 1. 



Remarks 7.1. 1. Notice here that we do not deal with the usual charge density and 
current density, since these ones are taken within the gyro- coordinates. 



2. We mention that we could have considered the more general case: 



g,{t,x,v)= / p^{t,x)l (^t .^^u{de)%^=o 



(7.14) 



M 



where (M, 0, i^) is a probability space which allows to model more realistic plasmas 
than "cold plasmas" and covers many interesting physical data, like multi-sheet elec- 
trons or water-bags data (we refer for instance to ^ and references therein). We will 
not do so for the sake of readability but we could deal with it with exactly the same 
analytic framework: the analogues of Theorems \2.1\ and \2.2\ identically hold. We get 
in the end the system: 
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dtvf + V±{E^vp + vfd^i {vp = x|| ) 

E^ = V^Al'{Jp^du-Jp^dx^du) ^ ' 

f p® {t, x)dx±du = 1. 

As before, the equations are coupled through x± and here also through the new pa- 
rameter Q. 

3. Actually, the choice: 

9e{t, X, V) = Pe{t, (7.16) 

leads to an ill-posed system. Indeed, we have to solve in this case equations of the 
form = fe,_L(t, X — v^) where v^^_\_ is the unknown. We can not say if this relation 
is invertible, even locally. 

Acknowledgments. This work originated from discussions with Maxime Hauray; I would 
like to thank him for that. 
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